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I.  INTRODUCTION 


A  muzzle  brake.  Figure  1,  is  used  to  reduce  the  recoil  of  the  gun 
thus  allowing  a  lighter  gun  cirriage.  The  impinging  gun-exhaust  plume 
exerts  counterforces  on  the  baffle  surfaces  thus  reducing  the  recoil 
impulse.  The  cowls  attach  the  baffles  to  the  gun  but  may  also  affect 
the  momentum  imparted  to  the  muzzle  brake.  The  objective  of  this  work 
is  to  calculate  the  muzzle-brake  momentum  index  by  a  simple  method. 

The  momentum  index  value  is  a  measure  of  the  counter  impulse  given  the 
brake  compared  to  the  gas  ejection  impulse  and  thus  is  a  measure  of  the 
brake  efficiency  or  effectiveness. 

This  work  stems  from  the  muzzle-brake  investigations  of  Oswatitsch1 ’2 
and  Smith3*1*,  Oswatitsch  shows  that  the  exhaust  flow  impinging  upon  a 
baffle  may  be  approximated  by  a  quasi-steady  jet.  Smith3’4  also  makes 
this  assumption.  Schmidt  and  Shear2  have  confirmed  this  hypothesis  by 
utilizing  optical  techniques.  Figure  2  displays  the  structure  of  the 
flow  from  the  muzzle.  Immediately  after  the  projectile  leaves  the 
muzzle,  the  gases  are  freed  to  expand.  The  free-air  blast  forms  and 
travels  rapidly  with  the  contact  surface  and  Mach  disc  following.  In 
the  supersonic  core,  the  gases  expand  from  the  high  pressures  at  the 
muzzle  to  very  low  pressures  near  the  barrel  shock.  Even  through  the 
Mach  disc  or  blast  front  is  moving,  the  region  in  the  supersonic  core 
approximates  the  flow  in  a  steady  jet  with  the  same  muzzle-exit  condi¬ 
tions.  With  the  baffle  placed  one  to  a  few  calibers  away  from  the 


1.  K.  Oswatitsch,  "Intermediate  Ballistics ,"  DVL  R  358,  Deutschen 
Versuchsanstalt  fur  Luft-unu  Baumfahrt,  Aachen ,  Germany, 

June  1964,  (UNCLASSIFIED).  AD  473  249. 

2.  K.  Oswatitsch,  " Flow  Research  to  Improve  the  Efficiency  of  Muzzle 
Brakes ,  Part  I-III, "  in  Muzzle  Brakes ,  Volume  II:  Theory, 

E.  Hammer  (ed.),  Franklin  Institute,  Philadelphia,  FA,  1949. 

3.  F.  Smith,  "Model  Experiments  on  Muzzle  Brakes,  "  RARDE  R  2/66, 

Royal  Armament  Research  and  Development  Establishment,  Font 
Halstead,  U.K.,  1966.  (UNCLASSIFIED) .  AD  48?  121. 

4.  F.  Smith,  "Model  Experiments  on  Muzzle  Brakes,  P;ct  III:  Measurement 
of  Pressure  Distribution,  "  RARDE  R  3/68,  Royal  Am  Trent  Research 
and  Development  Establishment,  Fort  Halstead,  U.K. ,  1868. 
(UNCLASSIFIED) .  AD  845  519. 

5.  E.  M.  Schmidt  and  D.  D.  Shear,  "Optical  Measurement  of  liuazie 
Blast,"  AIM  J . ,  Vol.  13,  No.  8,  August  1975,  pp.  1084-1'.  91. 
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muzzle,  the  baffle  would  be  in  the  supersonic  core  for  most  of  the  gun 
exhaust  phase.  Thus,  the  core  flow  impinging  uj on  the  baffles  primarily 
determines  the  momentum  index. 

In  the  present  work,  we  first  apply  a  simple  theory  for  predicting 
pressures  on  the  baffles  given  the  impinging  flow.  W  th  Schmidt,  et  al6, 
and  Smith3’4,  we  utilize  Newtonian  theory.  Next,  we  construct  a  simple 
model  for  the  flow  in  the  muzzle-jet  plume  assuming  a  supersonic  source 
flow.  With  these  assumptions,  we  then  vary  the  source  location  so  that 
the  theoretical  overpressure  curve  can  be  fitted  to  the  experimental 
data  of  Schmidt,  et  al6.  With  the  optimum  source  location  established, 
the  expression  for  the  pressure  can  be  integrated  and  a  simple  closed- 
form  expression  for  the  momentum- index  of  the  brake  is  obtained. 

The  resulting  model  is  applicable  to  most  gun  systems  although 
the  present  study  does  not  take  into  account  how  the  efficiency  varies 
with  the  amount  of  enclosure  associated  with  the  cowls.  However, 
Smith’s3’4  experiments  show  the  momentum- index  changes  little  for  small 
amounts  of  enclosure.  Figure  3  shows  that  when  5,  the  angle  that  the 
disc  is  inclined  from  the  radial  plane,  is  equal  to  zero,  (flat  discs), 
the  momentum  index  increases  monotonically  and  a^ost  linearly  as  the 
enclosure  wlue  for  the  cowls  is  increased.  When  che  cowls  enclose 
one-half  of  the  surrounding  space,  the  momentum- index  increase  is  about 
10%.  When  the  angle  6  is  greater  than  20° >  The  momentum  index  values 
go  both  below  and  above  the  value  for  the  fully  open  brake.  This  maxi¬ 
mum  excursion  when  surrounding  space  is  less  than  half  closed  by  cowls 
is  a  little  more  than  5%. 


If.  MODEL  FOR  THE  PRESSURE  AND  IMPINGING  FLOW  ONTO  THE  BAFFLE 

The  pressure  on  the  baffle  surface  is  calculated  by  Newtonian  flow 
theory.  Consider  that  the  flow  is  incident  on  the  baffle  surface  at 
the  angle  0  and  with  the  local  fluid  velocity  V.  As  fluid  strikes  the 
surface,  the  perpendicular  component  of  momentum  is  dissipated  but  the 
velocity  component  tangent  to  the  surface  is  unchanged.  The  resulting 
overpressure  equation  can  be  shown  to  be 

2  9 

p  -p  -  pV  sin~0  (1) 


6.  E.  M.  Schmidt ,  E.  J.  Gion ,  and  K.  S.  Pansier ,  "Measurement  of 
Muzzle  Blast  and  Its  Impingement  Upon  Surfaces , "  AIM  22th  Fluid 
and  Plasma  Dynamics  Conference ,  Williamsburg 3  VAt  July  1979 , 

AIM  Paper  79  -  1550. 
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Here  p  is  the  pressure  on  the  baffle  and  p  is  the  local  value  of  the 

pressure  in  the  flow  field.  Schmidt,  et  al6,  applied  this  theory  for 
gun-exhaust  flow  impinging  upon  a  disc-baffle  and  obtained  good  agree¬ 
ment  with  experiment.  They  approximated  the  supersonic-core  flow  onto 
the  baffle  as  a  steady-jet  method-of-characteristics  calculation 7 . 

Smith3  considers  "that  the  flow  which  would  normally  pass  through  the 
disc  is  simply  removed  from  the  flow  picture  with  a  corresponding 
removal  of  Its  thrust  component."  This  view,  for  a  straight  disc, 
is  equivalent  to  the  application  of  Newton's  flow  theory.  Nevertheless, 
for  conical  baffles  (6  >  0),  the  two  resulting  expressions  are  somewhat 
different. 

As  discussed  earlier,  the  flow  onto  the  brake  can  be  approximated 
as  a  steady  supersonic  jet-core  flow  with  muzzle-exit  conditions  for 
the  time  of  interest.  Figure  4  shows  a  steady  supersonic  core  flow 
generated  by  a  me^hod-of-characteristics  calculation7.  The  lines 
going  away  from  the  muzzle  are  streamlines.  The  other  set  of  lines  are 
isoproperty  lines  along  which  the  pressure,  temperature  and  Mach 
number  are  constant.  A  feature  of  this  flow  structure  permits  a 
simplifying  approximation.  A  short  distance  downstream  of  the  muzzle, 
the  streamlines  tend  to  straighten  out  although  still  possessing  some 
small  curvature.  Such  a  flow  is  reminiscent  of  the  simpler  supersonic 
source  flow, which  would  have  straight  streamlines  emanating  from  the 
center  with  the  isoproperty  lines  represented  by  concentric  circles. 

The  source  flow  is  defined  if  the  position  of  the  origin  and  the  gas- 
dynamic  properties  at  a  given  radius  are  known. 

Using  Newtonian  flow  theory  and  the  source  flow  approximation,  we 
can  derive  a  simple  expression  for  the  pressure  on  the  surface  of  a 
baffle.  We  can  start  by  restating  that  the  flow  field  of  a  source 
flow  can  be  determined  if  we  know  the  position  of  the  source  point 
and  the  gasdynamic  properties  at  a  given  radius.  We  will  use  the 
centerline  properties  to  find  the  flow  characteristics  elsewhere  m 
the  field.  We  use  a  method  of  characteristics  (MOCj  scheme  to  generate 
the  Mach  number  distribution  along  the  boreline  once  and  for  all.  Tnis 
distribution  approximates  sufficiently  well  the  distributions  for  a 
wide  range  of  muzzle-exit  conditions  of  interest.  Figure  5  shows  the 
Mach  number  distribution  and  other  centerline  properties  for  the  jet 
core.  Here  the  subscripts  j  and  c  denote  the  muzzle  exit  conditions 
and  the  local  boreline  conditions  respectively.  The  symbols  M  and  y 
denote  the  Mach  number  and  the  specific  heat  ratio  respectively. 


7.  A.  R.  Vick 3  E.  H.  Andrews ,  J.  S.  Pennard ,  and  C.  B.  Craidon , 

"Comparison  of  Experimental  Free-Jet  Boundaries  with  Theoretical 
Results  Obtained  with  the  Method  of  Characteristics ,  "  NASA 
Technical  Note  E-2327,  June  1964  (NTIS  N64-23032). 


The  distance  along  the  centerline  is  expressed  in  calibers.  In  fact, 
the  distances  in  this  paper  will  always  be  expressed  in  calibers. 

With  the  Mach  number  distribution  and  the  muzzle  exit  conditions  known, 
we  can  calculate  the  gasdynamic  properties  at  any  selected  point  on 
the  centerline. 


We  then  wish  to  know  the  local  gasdynamic  quantities,  such  as 
pressure,  velocity,  etc.,  located  on  a  baffle  surface.  We  first 
construct  two  spherical  surfaces  with  the  assumed  source  point  as 
their  centers.  One  surface  intersects  the  impingement  point  and  the 
other  surface  intersects  the  projection  of  the  baffle's  surface  upon 
the  boreline.  The  control  surface  is  defined  by  the  intersection  of 
these  surfaces  with  a  differential  cone  having  its  vertex  at  the 
source  point  and  its  axis  intersecting  the  impingement  point.  The 
local  gasdynamic  quantities  at  the  surface  intersecting  the  impingement 
point  can  be  found  by  the  momentum  theorem  since  the  local  gasdynamic 
properties  are  known  on  the  other  spherical  surface  from  an  MOC  calcu¬ 
lation.  Utilizing  an  approximation  to  the  Newtonian  pressure. 

Equation  (1),  an  expression  can  then  be  obtained  for  the  pressure  on 
the  baffle  surface. 


The  foregoing  brief  description  of  the  derivation  is  now  described 
in  more  detail.  Figure  6  shows  the  upper  cross-section  of  a  baffle 
and  the  end  of  the  gun  barrel.  The  origin  of  the  source  flow  is  along 
the  boreline  at  a  certain  distance  from  the  muzzle  as  shown  in  Figure  6. 
The  angle  6  is  the  angle  that  the  disc  is  inclined  from  the  radial  plane, 
<j>  is  the  angle  between  the  boreline  and  the  flow  direction  at  the  baffle 
surface  and  0  is  the  angle  of  incidence  of  the  impinging  flow.  With 
the  source  flow  assumption,  the  circle  defined  by  r  is  a  property 
contour.  In  order  to  find  the  pressure  at  the  given  impingement  point, 
we  can  construct  a  control  volume  from  a  differential  area  at  the 
impingement  point  to  a  differential  surface  at  r  as  shown  in  Figure 
This  control  volume  then  comprises  part  of  a  differential  cone  volume 
extending  from  the  source  point  to  the  impingement  point  on  the  baffle. 
The  following  relation  is  obtained  for  the  momentum  flux  tensor  in  a 
control  volume  with  steady  flow: 


pV.  V. 
1  3 


n . )  dS 
3 


0 


(2) 


Here,  standard  indicial  notation  is  used  for  sectors,  and  n.  is 
the  ith  component  of  the  unit  vector  normal  to  the  surface.1  For 
the  component  along  the  flow  direction.  Equation  (2)  would  yield 


(Pc  +  PCVC2)  dS5  =  Cp  +  pV2j  dS5  +  pm  (dS6  -  dS5)  (3) 

Here  the  subscript  c  denotes  the  flow  conditions  along  the  centerline 
at  Zy  and  pm  denotes  some  value  of  p  between  r  and  Tj.  The  last  term 

on  the  RHS  of  Equation  (3)  results  from  performing  the  integration  over 
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it 


f  i 

'$  ! 


the  sides  of  the  control  surface.  Since  dS  is  proportional  to  r  and 
by  the  law  of  sines,  r^/sin  a  =  r  /sin  0: 


dSr 


sin26 


dS^  sin2a 

Equation  (3)  becomes,  using  the  above  equation  and  neglecting  the 
relatively  small  RHS  second  term: 


,,2  ,  ,,  2,  sin  6 

p  +  pV  =  (pc  +  Pc  Vc  )  — y 

sm  a 

Now  in  the  flow  of  interest,  we  could  approximate  Equation  (1)  by 

Ps  =  (P  +  PV2)  sin20 


(4) 


(5) 


For  baffle  distances  greater  than  1  caliber  from  the  muzzle,  the  error 
in  the  approximation  to  Equation  (1)  is  small  for  the  baffle  configu¬ 
rations  of  interest.  Substituting  Equation  (4)  into  Equation  (5),  we 
obtain  in  terms  of  <p  and  6, 


Ps  =  (Pc  +  Pcvc2)  cos4  C4>-6)  /cos2S 


(6) 


To  transform  Equation  (6)  to  an  easily  calculated  form,  divide  through 
by  pM.  Multiplying  and  dividing  through  by  Pj ,  we  obtain 


!s  .Psfi 

p  p-  p 

r  CO  x  J  X  CO 


(7) 


The  value  of  Pj/P^  is  estimated  by  interior  ballistic  theory.  In  terms 

of  the  Mach  number  of  the  flow.  Equation  (6)  becomes,  using  Equation  (7) 
together  with  the  ideal  gas  law  and  the  isentropic  relations: 


(1  *  iiIm.V-1 

P  p  2  i 


where  M^  is  the  jet-exit  Mach  number  and 


cos  (<f>-6) 

2. 

cos  6 


(8) 
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Y-l 


(9) 


1  +  yM 


F  h 


Y-l  n 

(1  +  -V-  M  ) 

v  2  c  1 


The  value  of  F  as  a  function  of  distance  m  calibers  along  the 
boreline  is  shown  in  Figure  8.  These  values  of  F  were  calculated  with 
a  method-of-characteristics  numerical  scheme.  It  is  found  that  F  has 
negligible  dependence  on  expected  values  of  the  muzzle-exit  conditions 
,  y  and  p./p^,  for  baffle  positions  of  interest.  For  this  analysis, 

F  will  be  assumed  to  be  independent  of  these  muzzle-exit  conditions. 

The  thickness  of  the  line  comprises  the  differences  observed  for 
Y  =  1.2  and  y  =  1.3.  These  values  of  y  should  cover  the  expected 
range  of  y  for  propellant  gases.  Thus  F  is  applicable  to  almost  all 
gun  systems. 

Utilizing  Schmidt  et  al^,  experimental  data,  the  source  point 
position  may  be  varied  to  obtain  optimum  agreement  with  Equation  (8). 
The  expression  for  r  as  a  function  of  z  is  the  following: 

r  =  O.lOlz2  +  0.528z  +  0.271  +  tan  <$  (10) 


Here,  is  the  projectile  hole  diameter  expressed  in  calibers. 

Figure  9a  gives  the  overpressure  distribution  on  a  disc  that  is  located 
1  caliber  from  the  muzzle.  The  experiment  used  a  disc  with  pressure 
gages  inserted  at  half-caliber  intervals  along  the  radius  of  the  disc. 
Comparisons  between  the  experiment  and  the  curve  fit  show  good  agreement. 
Figure  9b  shows  the  comparison  at  1.5  calibers.  Figure  9c  giv>s  the 
comparison  for  the  baffle  positioned  at  4  calibers.  Here  the 
agreement  is  not  as  good.  It  is  impossible  to  move  the  source  to  a 
position  to  give  a  better  fit.  The  calculated  pressures  are  lower 
at  the  smaller  radii,  which  agrees  with  the  observations  of  Schmidt 
et  al6. 


Although  the  model  has  been  applied  to  a  20mm  gun,  the  universality 
of  the  plume-flow  characteristics  permits  the  application  of  Equation  (8) 
to  most  weapons.  The  model  would  appear  to  be  less  applicable  as  one 
encounters  higher  jet- exit  Mach  numbers.  From  the  Mach  contour  lines 
and  streamlines  generated  by  the  MOC  calculations,  the  flow  patterns 
indicate  that  the  equivalent  source  would  need  to  be  moved  farther  from 
the  muzzle  as  the  jet-exit  Mach  number  increased.  Compared  to  the 
sonic  jet,  this  would  lower  the  pressures  and  also  increase  the  fall- 
off  rate  of  the  pressures  with  increasing  radius. 
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III.  DEVELOPMENT  OF  MUZZLE-BRAKE  MOMENTUM- INDEX  EXPRESSION 


As  discussed  in  the  Appendix,  we  have  chosen  the  momentum  index 

to  be 


n 


(11) 


Here  I  is  the  recoil  impulse  of  the  gun-barrel  with  a  bare  muzzle,  I/ 
is  the  recoil  impulse  of  the  gun-barrel  with  the  muzzle  brake  and  G  is 
the  impulse  given  the  bare-muzzle  gun  by  the  propellant  gases.  The  gas 
impulse  is  essentially  the  time  integral  of  the  momentum  flux  passing 
the  muzzle  face  of  the  gun.  During  the  major  portion  of  the  gun¬ 
emptying  phase,  the  ratio  of  the  numerator  to  the  denominator  is  approxi¬ 
mately  a  constant,  as  shown  in  Appendix  A.  Thus,  Equation  (11)  is  well 
approximated  by  replacement  of  the  time  integrals  by  their  corresponding 
integrands. 


n 


4 


(p  -p  )  n  dS 
rs  z 


fp.+p .V.  )A. 
'yJ  3  3  3 


(12) 


Here  dS  is  a  differential  surface  on  the  muzzle  brake  and  n  is  the  z 
component  of  the  unit  vector  normal  to  the  surface.  The  denominator 
of  Equation  (12)  is  the  thrust,  T,  generated  by  the  propellant  gas  on 
the  weapon  and  the  numerator  is  simply  T-T'  where  T/  is  the  thrust  of 
the  propellant  gases  with  the  muzzle  brake  in  place.  Equation  (12) 
is  then 

T1 

•1  =  1-4  . 


If  l'  =  0,  then  all  the  exhaust  thrust  is  canceled  and  n  =  1.  The 
momentum  index  thus  gives  a  measure  of  the  effectiveness  or  efficiency 
of  the  brake.  Larger  values  of  n  indicate  a  more  effective  brake. 

From  Figure  6,  we  can  show  that  dS  in  terms  of  4>  and  r  is 

dS  =  (2rr2  cos26  sin<)>  d<j>)/cos3(<{>-6)  (13) 


Substituting  Equations  (8)  and  (,13)  into  Equation  (12)  with  the  value 
of  p  assumed  small,  we  obtain 

A  00  9 


n 


5  ^sinct 


2H  cos6  /  sincj>cos  (<{>-6)d<j> 


(14) 


where 


H 


4r2  (1  +  ^4  M.2)Y_1  F 
2  3 

1  +  yM.2 
3 


(15) 
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Integration  of  Equation  (14) 
n  =  H  cos<5  | 


yields  the  momentum  index  expression 
2  2 

cosd  (sin  d>  -  sin  d>  )  + 


sinS  [((j>m  -  cos<J>m  sin«J>m)  -  (^  -  cos^  sin<j>.  (16) 


Here  the  subscript  "m"  refers  to  the  maximum  value  of  rhe  angle  and  "n" 
refers  to  the  minimum  value  of  the  angle  subtended  by  the  muzzle  brake. 
Of  course,  the  total  force  on  the  baffle  can  be  obtained  by  multiplying 

2 

n  by  p.  (1  +  yMj  )Aj . 

In  keeping  with  the  spirit  of  the  preceding  approximations  and 
in  view  of  the  accuracy  expected  from  this  approach,  a  sufficiently 
good  approximation  to  the  calculated  H  is  achieved  using  the  following 
equation 


H  =  0.22  (z  +  R  tan6)  +  1.07  (1") 

The  value  of  H  is  quite  sensitive  to  the  prescription  for  the  source 
point  placement.  If  the  source  point  is  moved  toward  the  fixed  baffle, 
the  value  of  H  is  decreased.  Nevertheless,  the  multiplicand  of  H  in 
Equation  (16)  would  be  increased , thus  partly  offsetting  a  possible 
error  made  in  source  placement. 

To  illustrate  how  to  use  Equation  (16)  to  calculate  the  momentum 
index,  a  program  in  BASIC  is  presented  in  Appendix  B. 

It  should  be  noted  that  F  can  now  be  expressed  in  terms  of  Z,  6 
and  the  muzzle  exit  conditions.  From  Equation  (15),  we  have 


(i  +  ym.2)  h 

.......  — .  — — «L.„.  _  _ 

4r2  (1  +  M.2)  Y_1 


(18) 


The  value  r  is  expressed  in  terms  of  z  and  6  according  to  Equation  (10) 
and  H  is  found  from  Equation  (17).  Thus  the  expression  for  the 
pressure.  Equation  (8),  can  now  be  found  without  resorting  to  the  use 
of  graphs  or  extensive  calculations. 
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IV.  DISCUSSION 


It  is  of  interest  to  see  how  far  this  analysis  may  be  extended. 

For  a  variation  of  y  from  1.2  to  1.3,  the  value  of  H  varies  1%  for 
sonic  jet-exit  gas  velocities.  This  result  is  obtained  using  a  method- 
of-characteristics  calculation  with  Equation  (15).  Furthermore,  the 
Mach  contours  and  lines  of  flow  inclination  vary  little  in  this  range 
of  y.  Thus,  we  should  be  able  to  use  Equation  (16)  for  1.2  <  y  <  1.3. 
The  value  of  H  varies  more  with  the  parameter  M..  The  value  of  H 


increases  by  a  factor  of  1.35  from  =  1  to 


1.8  and  the  optimum 


source-point  location  for  agreement  with  experiment  probably  also  changes. 
To  see  how  the  momentum  index  varies  with  M.,  we  can  utilize  the  approach 
of  Schmidt6  which  describes  the  flow  by  HOC** calculations.  Figure  10 
shows  a  comparison  of  n  for  different  values  of  Mj  at  z  =  1.5.  The 
efficiency  is  found  to  vary  little  for  the  two  values  of  M.  for  discs  of 


different  radii.  Thus,  the  model  can  at  least  be  extended  to  the  M. 

range  1  to  1.8  for  baffles  in  the  vicinity  of  z  =  1.5.  Although  the 
insensitivity  of  H  to  undoubtedly  occurs  at  other  values  of  z,  there 

may  be  less  agreement  than  was  displayed  in  Figure  10.  For  most  cases, 
the  value  of  probably  declines  to  near  1  early  in  the  exhaust  phase, 

so  that  the  properly  weighted  value  for  in  Equation  (12)  would  be 

close  to  1.  Taking  into  account  this  canceling  behavior  and  the 
weighted  val  it  for  ,  the  error  should  be  small  in  assuming  a  sonic 

jet  for  most  cases. 


The  model  can  also  be  compared  to  experiment  for  varying  baffle 
angles.  Figure  11  is  a  comparison  between  Smith's3  experiments  and 
the  theory  for  baffles  placed  at  z  =  1.5.  The  efficiency  increases 
with  the  baffle  angle  partly  for  the  reason  that  the  baffle  is  inter¬ 
cepting  more  of  the  flow.  Smith's  experiment  was  run  for  y  =  1.67 
but  Smith  indicates  that  the  results  should  be  near  the  propellant-gas 
results  where  y  =  1,25. 


It  appears  that  the  source-flow  model  can  give  approximate  pre¬ 
dictions  of  the  momentum  index  for  different  design  parameters  of 
interest.  Figure  12  explores  the  influence  of  the  baffle  angle,  6, 
on  the  momentum  index  for  different  values  of  baffle  position.  The 
outer  and  inner  baffle  radii  are  held  constant.  The  momentum  index 
increases  little  for  larger  values  of  6  at  z  =  4.  In  fact,  the  momentum 
index  decreases  slightly  from  6  =  30°  to  6  =  40°. 

It  is  of  interest  to  find  the  distance  for  a  certain  diameter  disc 
that  would  produce  the  highest  momentum  index.  This  might  not  be  the 
best  distance  for  placement  since  the  brake  mass  and  strength  also 
have  to  be  taken  into  consideration.  Figure  13  shows  how 
n  varies  with  z  for  different  radii  of  discs  using  the  present  model. 


\ 
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The  model  results  are  compared  with  the  experimental  results  of  Smith. 
For  the  limited  amount  of  data  available,  the  agreement  between  experi¬ 
ment  and  theory  is  good.  Both  the  model  and  experiment  show  that  n 
decreases  as  z  decreases  for  small  values  of  z.  This  occurs  because  an 
increasing  amount  of  the  flow  is  going  through  the  projectile  hole. 

At  large  values  of  z,  less  of  the  flow  intercepts  the  brake  baffle  as 
z  increases.  The  maximum  value  of  n  occurs  at  larger  values  of  z  as 
the  disc  radius  increases.  The  effect  of  the  projectile  hole  loss 
becomes  apparent  at  larger  values  of  z  as  the  radius  of  the  disc  is 
increased.  The  maximum  values  for  the  experimental  curves  occur  at 
slightly  larger  values  of  z  than  the  model  predicts. 


V.  SUMMARY  AND  CONCLUSIONS 

A  simple  method  for  calculating  the  pressures  on  the  baffle  and 
the  impulse  given  the  muzzle  brake  during  firing  has  been  developed. 

The  method  utilizes  Newtonain  flow  theory  and  a  source  flow  whose 
origin  is  moved  with  baffle  distance  to  provide  optimum  agreement  with 
experimental  results.  Although  the  method  was  developed  using  one 
weapons  system,  the  model  can  be  applied  to  most  gun  systems.  Appli¬ 
cation  of  the  method  requires  knowledge  of  the  size  of  the  disc  and  pro¬ 
jectile  hole,  distance  of  the  projectile  hole  from  the  muzzle  and  the 
angle  6.  Additionally,  to  obtain  baffle  pressures,  knowledge  of  the 
muzzle-exit  conditions  is  required.  The  momentum  index  is  calculated 
using  a  simple  formula.  The  results  agree  well  with  an  independent 
experiment3' 4 „ 

Although  the  method  was  developed  specifically  for  disc  baffles 
and  sonic  muzzle  jets,  the  method  can  be  extended  to  more  complex  muzzle- 
brake  shapes  and  often  to  muzzle  flow  with  initial  supersonic  exit 
conditions.  Some  baffle  shapes  could  be  approximated  by  a  collection 
of  discs  and  sections  of  discs.  For  disc  baffles.  Smith3  indicates 
cowling  enclosing  the  surrounding  space  by  less  than  one-half  seems  to 
change  the  momentum  index  by  only  a  small  amount.  Nevertheless, 
comparing  the  present  method's  calculations  with  results8  for  rectangular 
baffle  plates  using  top  and  bottom  plates,  we  find  that  the  present 
method's  momentiun  index  values  are  lower.  The  average  values  obtained 
by  the  present  method  need  to  be  multiplied  by  1.3  to  agree  with  the 
average  experimental  values.  It  is  assumed  that  this  discrepancy  is 
caused  by  the  effect  of  the  top  and  bottom  plates. 
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1.2r-  Vo  =  VALUE  OF  7?  FOR  FULLY  OPEN  BRAKE 

Vf  =  VALUE  OF  rj  FOR  PARTIALLY  CLOSED  BRAKE 

f  =  FRACTION  OF  CIRCUMFERENCE  CLOSED  TO  FLOW 
-  S  =  0* 


Figure  3.  Effect  of  Partially  Closing  a  Brake 


Figure  4.  Streamlines  and  Isoproperty  Lines  for  Steady  Jet 
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Figure  5.  Some  Centerline  Property  Distributions  for  the  Jet  Core 


Figure  6.  Illustrating  Source  Flow  Impingement  onto  Baffle 


19 


Illustrating  Cone-Section  Control  Vclime  Located  next  to  the 
Baffle  Surface 
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Figure  8.  Jet  Property  Function  on  Borelir.c 
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A  ~  SCHMIDT,  GION  AND  PANSIER 
DISC  IS  1  CALIBER  FROM  MUZZLE 
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Figure  9a.  Comparison  of  Optimized  Model  with  Experiment  (z  =  1  caliber) 


Figure  9b.  Comparison  of  Optimized  Model  with  Experiment  (z  =  1.5  calibers) 
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A ~  SCHMIDT,  GION  AND  FANSLER 
DISC  IS  4  CALIBERS  FROM  MUZZLE 
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Figure  9c.  Comparison  of  Optimized  Model  with  Experiment  (z  =  4  calibers) 


Figure  10.  Momentum  Index  Comparison  for  Different  Muzzle-Exit  Mach 
Numbers 


Figure  11 . 
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Figure  12. 
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Momentum  Index  versus  Baffle  Angle;  Comparison  Between 
Theory  and  Experiment 


Momentum  Index  versus  Baffle  Angle  for  Different  Baffle 
Distances 
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APPENDIX  A.  MUZZLE  BRAKE  EFFECTIVENESS  DEFINITIONS 


The  muzzle-brake  effectiveness  parameter  gives  a  measure  of  the 
recoil  reduction  of  the  muzzle  brake.  There  are  various  definitions 
for  muzzle  brake  effectiveness.  Arguments  have  been  made  concerning 
the  best  definition  to  use.  In  this  Appendix  we  attempt  to  show  the 
advantages  of  our  definition. 

We  start  by  considering  the  momentum  given  the  gun  during  the 
firing  process.  We  neglect  friction  processes,  and  write  Newton's 
Second  Law  as 


3t/’Vi  d,  .  -/pn.dS  CA-D 

This  equation  is  in  indicial  notation.  Here  dx  is  a  volume  element, 
dS  is  an  element  of  the  surface  chosen  and  n^  is  a  component  of  the 

unit  vector  normal  to  and  pointing  out  from  the  control  surface.  This 
expression  can  be  manipulated  to 

M  +  i_  /pV.  dx  =  -  /( pn.  +  pV.V.n.  )  dS  (A-2) 

p  dt  p  3tJ  1  J  i  i  k  kJ 

Here  Mp  and  V  are  the  mass  and  velocity  of  the  projectile  respectively. 
Placing  the  control  surface  along  the  interior  surface  of  the  gun  and 
across  the  surface  of  a  bare  muz  :le  gun,  we  obtain'  for  times  before  shot 
ejection 


>bdSb  -/Pdssnz  =  /(pj 


+  p.  V.  )  dA.  +  ~ 
ii  i  at 


pV  dx  +  M  V  (A-3) 
z  p  dt  p  v  ' 


Here  the  subscript  j  denotes  exit  conditions,  V  is  the  z- component 

of  the  velocity,  the  subscript  b  refers  to  the  interior  base  of  the 
gun  barrel,  the  subscript  S  refers  to  the  sides  of  the  tube  and  dA^ 

refers  to  the  element  of  surface  at  the  gun  muzzle.  The  propellant 
gas  and  the  precursor  gas  are  included  in  the  element  dx  in  the  second 
term  on  the  RHS  of  Equation  (A-3).  Now  assume  that  the  gas  dynamic 
conditions  do  not  vary  radially.  The  total  recoil  impulse  can  be 
obtained  by  subtracting  pM  from  both  sides  of  Equation  (A-3)  and 
integrating  with  respect  to  time: 

+00 

r  +  03 

1  =  Aj  j  (Pj-Poo  +  pjvj2)  dt  +  mp  vp  +fftt  (pvdx)dt 


Here,  t  =  0  corresponds  to  the  time  at  shot  ejection.  The  last  term. 
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when  evaluated  at  t  =  +  »,  drops  out.  The  impulse  for  the  gun  with 
the  muzzle  brake  is 

=  I  -f  [jlvs-  Pj  nz  dS]  dt  (A-5) 

where  the  integrand  is  integrated  over  the  brake  surface,  S,  and  ps  is 
the  pressure  on  the  brake  surface. 

The  momentum  index  can  be  defined  as 


_  I-l' 
n  =  — 


(A-6) 


where  G  is  the  impulse  given  the  gun  by  the  propellant  gases 


G  = 


I  - 


m  V 
P  P 


(A-7) 


Thus,  the  momentum  index  becomes 


n 


dS]  dt 


G 


(A-S) 


The  flow  pattern  in  the  jet  varies  little  with  the  jet-exit  pressures 
of  interest.  Over  the  range  of  jet-exit  pressures  of  interest,  the 
integrand  part  of  the  numerator  should  be  proportional  to  the  muzzle- 
exit  pressure  and  the  integrand  in  the  denominator  is  also  proportional 
to  the  jet-exit  pressure.  Furthermore,  results  obtained  by  MOC 
calculations  and  application  of  Newtonian  theory  show  that  the  ratio 
of  the  integrands  in  Equation  (A-8)  is  only  weakly  dependent  upon  the 
initial  jet-exit  Mach  number.  The  momentum  index  can  then  be  approxi¬ 
mated  by 

n  = 


(p.  +  p.v.  )  A. 
*1  JJ  J 


The  exit  properties  would  be  weighted  so  that  they  represent  times 
shortly  after  shot  ejection  when  muzzle  pressures  were  comparable  to 
the  peak  muzzle  pressure.  The  flow  pattern  should  be  only  very  weakly 
dependent  on  charge  and  projectile  masses.  The  value  of  n  should  also 
be  very  weakly  dependent  on  the  initial  muzzle  exit  conditions.  Hence, 
n  should  depend  primarily  on  the  muzzle-brake  configuration. 
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The  experimental  determination  of  n  is  rather  straightforward. 

One  simply  determines  I  and  l'  with  an  impulse  measuring  machine. 

The  projectile  mass  and  velocity  of  the  projectile  can  also  be  measured 
with  some  accuracy.  The  value  of  G  is  easily  determined  by  liquation 
(A- 7 ) .  The  value  of  the  momentum  index  is  then  determined  by  substi¬ 
tuting  experimental  values  into  Equation  (A-6).  The  expression  for  n. 
Equation  (A-9) ,  also  lends  itself  well  to  theoretical  analysis. 

L.  L.  Pater8  has  defined  a  similar  momentum  index  except  that  G 
is  replaced  by  I  ,  "defined  as  that  portion  of  the  total  recoil 

impulse  which  occurs  after  projectile  ejection."  This  value  of  the 
impulse  would  be,  using  Equation  (A-4)  and  (A-7)  but  evaluating  at  t  =  0 
instead  of  at  t  =  +°°, 

I  =  G  -  A.  f  pVdz  -  A.  f  (p.  -  p  +  p.V.2)  dt. 
a  iy_L  J  J  00  J  J 

Here  L  is  the  length  of  the  barrel  in  calibers.  The  third  term  in¬ 
volves  the  flow  of  precursor  gases  out  of  the  barrel  and,  in  most 
cases,  would  be  small  compared  to  the  other  terms.  There  are  two 
difficulties  with  using  I  .  First,  one  must  make  an  assumption  about 

cl 

the  distributions  of  p  and  V  inside  the  barrel  at  the  time  of  shot 
ejection  since  they  are  not  known.  This  can  lead  to  some  experimental 
and  theoretical  uncertainties.  Second,  for  this  definition  of  momentum 
index  to  be  invariant  to  charge  weights,  projectile  weights  and  barrel 
length,  the  second  and  third  terms  must  be  directly  proportional  to  the 
corresponding  values  for  G.  It  is  not  at  all  obvious  to  the  author 
that  this  would  be  the  case  although  L.  L.  Pater  noted  that  in  his 
experiments,  his  momentum  index  for  a  given  muzzle  brake  varied  little 
with  the  different  gun  parameters  he  used. 

There  are  other  definitions  for  measures  of  muzzle  brake  effective¬ 
ness.  L.  L.  Pater8  discusses  these  in  some  detail  in  his  report.  He 
shows  that  they  all  seem  to  have  some  drawbacks  that  limit  their  use¬ 
fulness  . 


8.  L.  L.  Pater ,  "Muzzle  Brake  Parameter  Study  j"  NSWC/DL  TP- 3531 , 

U.S.  Naval  Surface  Weapons  Center 3  Dahlaren  Laboratory ,  Dahlqren , 
VA  22448 ,  October  1976 ,  ( UNCLASSIFIED )\ 


APPENDIX  B.  PROGRAM  FOR  CALCULATING  THE  PRESSURE  AND  THE 
MOMENTUM  INDEX 

Below  is  a  computer  program  coded  in  BASIC  that,  given  the  pro¬ 
jectile  hole  radius  and  outer  baffle  radius,  can  calculate  the  momentum 
index  either  as  a  function  of  baffle  distance  from  muzzle  or  the  baffle 
angle,  6.  It  can  also  calculate  the  pressure  field  on  a  baffle. 


1  i'o  t>Ll(  PR  OCR  Mi  In  f  I  HD  II  OH  El  ITU  II  INDEX  OR  PRESSURE  Oil  6ATFLE 
lie  PRINT  MO  '  .  i,  UhHT  TO  ChLCULhTE  THE  PRESSURE  '< " 

1  -  e  1 1  i  P  U  T  E  4 
1  So  IF  =  '  TtlF.II  660 

HO  FF1HT  "INHJ7  THE  PROJECTILE  HOLE  RADIUS  IN  CALIBERS" 

150  REIT  R 1  IS  SHhLLEST  RADIUS 
led  INPUT  R 1 

170  REM  RO  IS  f HE  MAXIMUM  RADIUS  OF  THE  BAFFLE 
180  PRINT  "INPUT  THE  BAFFLE  RADIUS  IN  CALIBERS" 

ISO  INPUT  RO 

200  PRINT  "DO  YOU  WhNT  TO  CALCULATE  MOMENTUM  INDEX  III  TERMS  OF  2" 

210  PRINT  "ANSWER  BY  PRESSING  EITHER  Y  OP  M" 

220  INPUT  h$ 

230  lr  Ai  = " 7  “  THEN  333 

240  PRINT  "INPUT  THE  DISTANCE > 2>  TO  BhFFLE  IN  CALIBERS" 

250  INPUT  2 

260  PRINT  "TO  START  ODER,  INPUT  THE  UALUE  100" 

270  PRINT  "INPUT  THE  ANGLE  IN  DEGREES" 

280  INPUT  D 

280  IF  D=>100  THEN  140 

300  1-1 

310  GO  TO  420 

320  REM  D  IS  THF  BAFFLE  DISH  ANGLE<DEGREES> 

3Ju  PRINT  "INPUT  THE  ANGLE  DELTA  IN  DEGREES" 

340  INPUT  D 

350  REM  2  IS  DISTANCE  FROM  MUZZLE  TO  RADIAL  PROJECTION  OF  R1  ONTO  AXIS 
360  PRINT  "TO  SThRT  OUER,  INPUT  THE  UALUE  100" 

370  PRINT  "INPUT  THE  DISTAHCE,Z,  TO  BAFFLE  IN  CALIBERS” 

380  INPUT  2 

350  IF  2=  > 1 00  THEN  140 
400  1=0 

410  REM  THE  ANGLE  "DELTA"  IN  TERMS  OF  RADIANS 
420  D1=PI*D/180 

430  REM  26  IS  A  RESULT  OF  AH  INTERMEDIATE  CALCULATION 
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PhECEDIMJ  page  BLANK' 


■NOT  FILMED 


4  40  :<f=0. 271  +  0. 528*2*0.  191*212 

450  R  r M  25  IS  A  DISTANCE  CI"EN  AS  "r"  IN  REPORT 

460  25  =  26+Pi*TAN*  D1 '• 

4 TO  PFH  mNGLE  IN  REPORT  CORRESPONDING  TO  PI 
4S0  r  !  =  A  fNkR  1  i  Z5-R1 t TAN<  D1 )  -■ ' 

490  REN  ANGLE  IN  RLPORT  CORRESPONDING  TO  R0 
500  PO  =  ATN'  RO  •  .15  -ROtTANCDl  V, 

510  C  1  =COS'  F 1  '• 

528  S1=SIN* F 1 j 
530  CO=CO3<F0  > 

540  S0  =  SIN(F8'* 

550  C2=C0$< D1 > 

560  S2~SIH''D1 ) 

570  REM  H  AS  DEFINED  IN  REPORT 
580  H=1.0?+Q.22t<Z+Rl*TAN<Dro 
590  T1=F1 -CltSl 
600  T0=F0-COtS0 

6I0  REI1  E  IS  THE  MOMENTUM  INDEX  AS  DEFINED  IN  REPORT 
620  E=HtC2*<X2t'’S0t2-S112>+S2t<T0-Tl>> 

63C  PRINT  D ,  Z  <  E 
640  IF  1=1  THEN  280 
658  GO  TO  380 

660  PRINT  "INPUT  THE  MUZZLE  EXIT  PRESSURE  IN  ATMOSPHERES” 

670  INPUT  PI 

680  PRINT  "INPUT  VUE  JET  EXIT  MACH  NUMBER" 

690  INPUT  Ml 

700  PRINT  "INPUT  DISTANCE  FROM  MUZZLE  TO  BAFFLE" 

710  INPUT  Z 

720  PRINT  "INPUT  THE  BAFFLE  ANGLE. DELTA,  IN  DEGREES" 

730  INPUT  D 

740  PRINT  "INPUT  THE  PROJECTILE  HOLE  RADIUS  IN  CALIBERS" 

750  INPUT  R1 

760  PPINT  "TO  START  OUER,  INPUT  THE  UALUE  100" 

770  PRINT  " INPUT  THE  RADIAL  POSITION  OF  THE  BAFFLE  IMPINGEMENT  POINT" 
780  INPUT  P 

79o  IF  R  =  Muu  1  HEN  1  Cm 

SGu  REM  THE  hNGLE  "DELTA'  IN  TERMS  OF  RADIANS 
SIC*  D  1  =  P  1 1 D  ISO 

S20  REI1  26  IS  h  RESULT  uF  AN  INTERMEDIATE  CALCULATION 

S30  26  =  0.271 '0.52612  +  0.  101*212 

840  REM  25  IS  h  DISTANCE  GIUEN  AS  V  IN  REPORT 

850  25  =  26  +  R 1 TTmN<  D1 » 

860  REM  ANGLE  "PHE"  IN  REPORT  CORRESPONDING  TO  R 
87o  PS=ATN*R  <  25-R*TAN(Dl > ) > 

880  REM  21  IS  AS  DEFINED  IN  REPORT 
890  :i  =  :-tR11TAN*Dl> 

900  REM  H  IS  AS  DEFINED  IN  REPORT 

910  H=0. 22*21+ 1 .07 

920  REM  F  IS  AS  DEFINED  IN  REPORT 

930  F=t . 251H- 4  -  25=25 

940  REM  INTERHEDIhTE  CALCULATION 

950  1 1  =  ‘  1+H1 1M1  ;'8*t5 

960  PEI!  EXPRESSION  FOR  PRESSURE  ON  SURFACE  OF  BAFFLE 
970  P=P1 1 1 1  !FlC08(R8-Dl  >  14 -'COS*  D1  >t2 
980  REM 
990  P2=P- 1 

1000  REM  THE  ANGLE  "PHE"  IN  DEGREES 
1010  P9=  160*P8-'PI 
1620  PEM 
1030  REM 

1040  PRINT  R,P9,F,P2 
1050  GO  TO  780 
1060  END 
107O  REM 
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LIST  OF  SYMBOLS 


•\  Area  of  gun  bore 

f  fraction  of  cylindrical  surface  defined  by  R  the  muzzle  plane, 

and  plane  intersecting  baffle  outer  radius  that  is  closed  by 
cowling 

F  dimensionless  factor  in  expression  for  baffle-pressure 

G  propellant-gas  impulse  given  gun  without  muzzle  brake 

11  dimensionless  factor  in  momentum  index  expression  that  is  dependent 
on  baffle  position 

I  recoil  impulse  without  baffle 

I  propellant-gas  impulse  defined  by  Pater 

B 

l'  recoil  impulse  with  baffle 

nip  mass  of  projectile 

M  local  Mach  number  in  flow 

M  centerline  Mach  number  at  z^  (refer  to  Figure  6) 

M.  Mach  number  at  jet  exit 

1 

Ps  pressure  on  baffle  surface 

Pc  jet  centerline  static  pressure  at  z^ 

p.  static  pressure  at  jet  exit 

p  local  pressure  in  jet 

Pm  ambient  pressure 

r  distance  from  source  point  to  the  intersection  point  of  the  baffle 
plane  with  the  boreline  expressed  in  calibers 

r^  distance  from  source  point  to  a  point  on  the  baffle  surface 
expressed  in  calibers 

R  radial  position  on  baffle  surface 

R  outer  radius  of  baffle  expressed  in  calibers 
o 
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LIST  OF  SYMBOLS  (Continued) 


R1  projectile  hole  radius  expressed  in  calibers 

5  surface  to  be  integrated 

t  time  v»ith  respect  to  projectile  uncorking 

T  thrust  given  to  gun  by  propellant  gases 

V  local  fluid  velocity 

jet  velocity  along  centerline  at 

Vp  velocity  of  projectile 

z  position  of  baffle  with  respecv  to  muzzle  expressed  in  calibers 

z^  distance  from  muzzle  to  the  intersection  point  of  baffle  plane 

expressed  in  calibers 

a  angle  that  the  baffle  surface  makes  with  the  boreline 

Y  specific  heat  ratio 

6  angle  that  the  disc  is  inclined  from  the  radial  plane 

p  local  fluid  density  in  jet 

local  fluid  density  along  centerline  at  z^ 

6  angle  of  incidence  for  gas  flow  onto  plate 

$  angle  between  boreline  direction  and  flow  direction  at  the 

baffle  surface 

n  momentum  index  defined  by  Equation  (11) 
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